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1 Introduction 



Quasi-exactly solvable (Q.E.S.) Schrodinger equations nowadays attract much attention. They 
are equations for which a finite number of solutions can be analytically determined. In a 
sense, they constitute the intermediate step between exactly solvable equations (associated 
with potentials such as the harmonic oscillator one, the Coulomb one,...) whose all analytic 
solutions can be obtained, and the analytically unsolvable ones requiring a numerical treatment. 
These Q.E.S. equations have been studied following two different points of view. 

- First, the number of solutions can be related [Q, |3|, ^ to the dimension of the irreducible 
representations of a Lie algebra, namely sl{2,R). In that case the corresponding Schrodinger 
Hamiltonian can be written as linear and quadratic combinations of the sl{2,R) generators, 
these ones preserving the finite-dimensional space of the solutions. Due to the preservation of 
this space, the associated time-independent differential Schrodinger equation actually reduces 
to an algebraic one. 

- Second, this reduction from differential to algebraic is maintained but the number of 
known solutions (essentially two fs], |6| or three 0]) is fixed at the start. Generally speaking 
this second point of view is not subtended by any Lie algebra except in the case of the two 
solutions Instead, it is relevant of supersymmetric quantum mechanics Q. 

The main purpose of this paper is to propose a general approach unifying these two points 
of view. Our approach, which is also algebraic and, for the moment, restricted to the real line 
] — oo, oo[, or half line ]0,oo[, is well-suited for all potentials leading to a solution given by an 
overall function modulated by a polynomial. To our knowledge, it reproduces all the known 
exact cases which can be expressed in terms of polynomials, and most, if not all, of Q.E.S. 
models. For example, for the exact cases, it reproduces solutions of the harmonic oscillator, 
the Morse and Poschl- Teller potentials, solutions of the Lame equation |10|, the even potentials 
of the class described recently by Bender and Wang in Ref. On the contrary, since the 

odd potentials of the latter class cannot be written in term of polynomials (the solution is 
written as a non-terminating confiuent hypergeometric function), they cannot be dealt with. 
As examples of Q.E.S. Schrodinger equations, we can quote all the potentials given by Turbiner 
in a variation on the Lame equation as described in as well as the Tkachuk potential 
of Ref. The list is non-exhaustive. 

The paper is organized as follows: In Section ^, we explain the formalism of our method in 
a simplifying limit regarding the expansion basis. In Section ^ we develop the most general 
algorithm based on a single basis function, and restricted to finite sums. In Section ^ we show 
the relevance of the method by investigating several well-known Q.E.S. equations, of both 
types quoted above. We show in Section |5| that it can be applied to yet unused potential by 
considering a modification of the sextic oscillator potential. Finally, we draw our conclusions 
in Section |6[ 
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2 A Q.E.S. general approach 



As stated in the Introduction, we concentrate on the one-dimensional and time-independent 
Schrodinger equation 



iI:n{x) = En^n{x), iV = 0,l,2,.. 



(1) 



for which a finite number [= n + 1) of solutions iPm{x) can be determined. The procedure we 
shall explain here is not the most general one, which can be found in Section ^. We proceed 
so because we want to give here more details relevant to understand the philosophy of our 
approach. In order to cover the Lie algebraic and non-Lie algebraic Q.E.S. equations we write 
the potential V{x) as 



V{x) 



Vmax 

E 

k=-2 



(2) 



ymax being a positive integer, while the coefficients and the function f{x) are, for the 
moment, arbitrary. Making the sum begin with the index k = —2 allows to take care of the 
centrifugal barrier in the radial equation of a 3-dimensional problem. We shall see in Eqs. (|ll|) 
and ( [I3D that the minimal value of k has anyway to be A; = —2 in order to factor out a common 
power of f{x). 

The eigenfunctions iPn{x) are assumed to be 



N 

/ ^ '-m, 
m=0 



Xx)Y 



(3) 



The function g{x) plays the role of a weight factor, [/(x)]'^^ is introduced due to eventual 
singularities in (|2|), and the quantities ch^^ are expansion coefficients on the basis h{x). We 
introduce the notation 



N 



with 

We then have 



^Pm{x)=g{x)[f{xt^[h{x)Y 



(4) 
(5) 



XNi^N — 1) 



i^mix) 

f'Hx) 

P{x) 



g"{x) ^ fix) h"{x) 



g{x) 



mim 



1) 



h'^x) 
'hHx) 



fix) 
2Xn 



m- 



h{x) 



f'{x)g'{x) 
f{x)g{x) 



2m 



g'{x)h'{x) 

g{x)h{x) 



•^mix), 



(6) 



where the prime stands for the derivative with respect to x. The key idea for ensuring the 
quasi-exact solvability of (||) is to stabilize the space of ipm{x) {m = 0, 1, A''), or in other 
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words to ask for the second member of (|6|) to be a linear combination of the 'ipm{x). Because 
of the different dependences of the terms of this second member (with respect to m, etc.), we 
ensure the quasi-exact solvabihty by requiring 



AI 



fix) 

fix) 



g'ix) 



h'{x) 



M 
1=0 

M 

-g{x)Y,9}[h{x)]\ 
1=0 

M 

Y.^}[Kx)t 



(7) 
(8) 

(9) 

(10) 



the different coefficients fi ,fi, g}, h\ having to be fixed according to the potentials we have to 
consider. The upper-right index indicates the derivative order of the corresponding expanded 
function, e.g., to f{x) correspond the expansion coefficients ff, while to f'{x) correspond the 
coefficients //. The upper summation index M might be different for f (x) , f (x) , g' (x) , h' (x) . 
For simplicity, we take the same value: It is simply the highest value of the three different 
indices, the added expansion coefficients being vanishing. For example, if the sum leading to 
f(x) in Eq. (|^) is over the range / S [0, P] with P < M, we can extend it to the range / G [0, M] 
with = for z G [P + 1, M]. Since h{x) is a function to be chosen from the beginning, /'(x) 
and h'{x) are known. One just has to consider M big enough in order to be able to determine 

gix). 

Taking care of (p[)-(p^), we can rewrite (^) as 



+ Vix) 



-m 



X gix) [fix)]'- [hix)r+^ + [ _ ^(i - i + l)h]flj^, 
+ 2Xr, f]gl_^ - 2m\N ^ h]fl^^,] g(x)[/(x)]^--i[/i(x)]™+^ 
+ [ - X^iX^ - 1) E f]f^ gix)[fix)]'--\hix)] 



+Y,yk9ix)[fix)] 



h{x)Y 



where we have used 



i=0 j=0 i=0 j=0 



:ii) 



(12) 
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and where, for simplicity, we have omitted all the summation indices^ They will anyway not 
be needed in the general procedure to be described in the next Section. For the moment, we 
have not yet considered the expansion (^). It is used in (|ll]) when extracting a common factor 

+V{X)\ ^m{x) = Y.^^Zl+^9{x)[f{x)]'--'[h{x)Y 



(13) 



with 



fO 

h 



€L+^ = -^^(^^ - 1) E f}f^j - -j-h + i)/.]/i,_,,+i/,° 

- E 9]gU-i,-i,f?J?, + E - J' - ^1 - ^2 + i)h]gl,^i,^i,+,flh 

j,h,h j,h,l2 

+ E ykfi^fi^---fif^^ji-i^-i2-...-i 

k,h,l2,...lk+i 



m[m 



1) ^ h^jh}^j_i^_i^^2fhfh 



l-k + l 



We thus have 



m+i 



which, according to (|l]), must be equal to 



m+i 



(14) 



(15) 



(16) 



Equating the respective second members of ( p^ ) and (|T6|), we are led to the final algebraic 
equations, valid for all possible values of j, giving rise to ipNix), as well as En- 



N 



0, 



(17) 



m=0 



where, let us recall it, hm,j is defined through (|T^. We have thus reformulated the problem (|l]) 
to the form (|l7|). The latter looks heavy to handle. However, it can be straightforwardly 
implemented into a symbolic computer program. 

Before turning to the general procedure, let us finally notice that the approach developed 
in this Section reduces to the one which can be found in if V-i = V-2 = 0, \n = and 
fix) = h{x). 

*The relation between the coefficients fi and ff of Eqs. (Q) and can be obtained with the help of Eqs. ( [lo| ) 
and @. This relation is // = J^k hkf°+i-k{l + 1 - fc)- 
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3 A more general approach 



In this Section, we briefly describe the most general algorithm we have developped starting from 
the approach of the previous Section. As far as we now, there is no extension of it if we restrict 
ourselves to one basis function h{x). A generalization to more basis functions /i2(x), ...) 

is possible. In particular, some of the potentials we have studied may be treated using this 
generalization. However, they may as well be treated using the algorithm we present below. 
For this reason, the study of the generalized algorithm based on more basis functions is left for 
a future work. 

In the previous Section, we have described in details the philosophy of our approach to find 
analytical solutions to a given Schrodinger equation. Those analytical solutions are possible if 
we can factor out appropriate powers of h{x) and f{x), thus closing the system, see Eq. (|l3|). 
This later equation was found using the expansion of the potential on a basis f{x), see Eq. (|2|), 
together with the assumptions (|^)~(|lO|). We show now how to relax these constraints in order 
to treat the maximum number of different potentials. 

The simplest extension is to allow negative indices for the sums in Eqs. This 
brings no difficulties, and this is why we did not specified the range of the indices in the 
previous Section: 

oo oo oo oo 

^ fli ^ 6j= ^ bja,^j = Y,Y,bjai^j, (18) 

i=— oo j=—oo j=— ooj=— oo i j 

where, from now on, an index without range means a running on Z. Taking care of negative 
indices is required for example when using a basis such as h{x) = exp(— x). One can understand 
that to obtain a solution starting like cosh(x) may require negative powers of h{x): 2cosh(x) = 
h{x) + h{x)~^. For the same reason, the expansion of (^) is extended to start from negative 
orders. 

To allow more flexibility in the choice of the potential V{x) in (^), it is also necessary to 
make the sum begining to a lower negative index than k = —2. For this reason, we make it 
begin at k = —Vmin, with the condition Vmm > 2 this restriction being due to the fact that at 
least a factor of f[x)^^~'^ must be factored out in order to close the system. 

In studying the way the different powers of h{x) enter into Eq. (^, we can also convince 
ourselves that integer powers of the square root of the right-hand-side of Eq. ( [lO[) can be chosen. 
The analysis can be done using an arbitray power a: h' [x) = {Yl'^-oo ^H^(^)]'} • Plugging 
this expansion into (|^), it is straighforward to see that the condition 2a = n G N is required. 
Moreover, using ( [T8| ) and properly redefining the expansion coefficients, we can always write 
this h'{x) as the product of a series of the form of ( |lO|) with a series of the form of the square 
root of (0). 

Using a square root in h'{x) also requires to generalize the term g{x) in order, for example, 
to treat properly the term g'{x)h'{x) in Eq. (^). The generalization consists in using the 
exponential of a /i(x)-series, instead of the differential form given in or, equivalently, to 
multiply the latter by h'{x) and properly redefine the expansion coefficients. 

To summarize, the different expansions replacing Eqs. (^), and (p[)~([lO|) are: 

N 

i,N{x) = Yl (19) 
m=-N 
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Vmax 



V{x) = Yl ^kifixt, (20) 

A;=— Vmin 

fix) = (21) 

g{x) = exp|-^5/[M^)]'-5iogln[M^)]|, (22) 

h'{x) = |^/i,MM^)]j|E^H^(^)]'} ' (23) 

where we have added a term g^^^ln[h(x)] in Eq. ([2^). This is equivalent to multiplying Eq. 
by /i(x)~^'°s, thus generalizing again the approach of the previous Section, and increasing its 
flexibility. We recover the algorithm of the previous Section by taking = 0,hl = 5ifi, 

9\og ~ 0' 9l ~ Yln^ndi-n+ii^ — n + 1). We also remark that not all the unknowns are 
independent: For example, a problem with = l,gl^g = implies the sum, see Eq. (|l9|) , 

ci^lh{x)~^ + Cq^^ + cl^h{x). This is equivalent to h{x)~^[ci^l + '^'o^h{x) + c^i^h{x)'^]. A similar 

result is achieved with = 2,^J^g = l,c^2 — ^-i — ^^'^'o^ — ^^-h'^^i^ = '^o^^'^2^^ — '^^i^- 

The last possible generalization we quote is to modify again the expansion of the potential. 
The function f{x) in Eq. ( pO| ) can be modulated by any linear combination of the basis function 
h{x). Indeed, with new coefficients vi^k of two indices refering to the power k of f{x) and to 
the power I of h{x), the potential 

ymax 

V{x)= V,[f{xtY%'^iHx)f^ (24) 

k=—Vmin I 

when multiplied by iprn{x) of Eq. (^) just leads to the factorization of a power of h{x). This is 
seen as follows: 

V{x)'4)m{x) = 
ymin+Vmax 

E Vk-Vm^nfixf Yl %k-Vmin[Hx)fg{x)[h{x)rf{x)^'^-^'-'- 
k=0 I 
ymin+Vmax 

= Y ^k-Vmin Y 'fX^T Y %k-Vmin[Hx)f9{x) [/l(x)]-/(x)^--^-- 
A:=0 n I 

ymin+ymax 

= Y ^fc-ymin5^'^°[/i(x)]'+™5(x)/(x)^--^--, (25) 
A:=0 / 

where we have used — compare with the last line of Eq. (IT^ ) — 
and where 
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I 



E'^°[^(^)]'' (27) 



the last equality defining the coefficient ^ff . 

In Eq. (26), the coefficients ^f^ are defined recursivelly: 

f{xf = fix) fix)'-' =^ 'ff = '-'f?-nfl (28) 

n 

with the starting coefficient ^ff = 5i^q. In defining the different coefficients, we have taken 
the convention that the upper-left index refers to the power of the given function, while, as in 
the previous Section, the upper-right index indicates the order of its derivative. (We do not, 
however, specify the upper-left index when it is equal to 1, e.g., V/ = //•) 

With the generalization given in Eqs. (21)-(^^, and with the definitions of the coefficients 
in (p7|)-(p8D, we are now prepared to give the full algorithm solving Eq. ([l|). For this, we 
introduce also the following notation for the different terms entering Eq. (|6|): 

(29) 
(30) 
(31) 
(32) 
(33) 

g'ix) = -gix)\^c,}^^il + l)[hix)]'+gljhix)Y'{x) 

^ -gix)h'ix)Y,g}[hix)]\ (34) 

¥ = ^gld-m^ (35) 

m 

[g'ix)f = [gix)]'Y.[hix)fY.'^l'gl^^gix)'Y.'gl[hix)]\ (36) 

/ m I 

f'ix)h'ix) = ^[Mx)]' = J^AHM^)]', (37) 

/ m I 

g'ix)h'ix) = -5(x)J][/i(x)]'E2;^U-™^ (38) 

I m I 

f'ix)g'ix) = -gix)[h'ix)fY.ihix)fJ29Uim = -gix)[h'ix)fY.M[hix)f 



'hi 


— / J "'m,"'l-mi 
m 




[h'ix)f 


= T.^Hx)fY'hlhl^^ 

I rn 


Y,'h\[hix)t 
I 


fix) 


= j;/,°+i(/ + l)Kx)]'/i'(x 


)^E/'[^(^)]'^'(^) 




/ J J mJ I— ml 
m 




[f'ix)f 


= Y.\Kx)]'Y."hl'flm^ 


E'//[Mx)]', 
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h"{x) 



fix) 



9"{x) 



= -g{x) Y,[h{x)f J2 fa'jhl^ ^ -9{x) fg][h{x)t 

I m I 

= \^ /,! 

m 

/ m 

I 

= Y.'^KX)]' Yl ^hlnll+l-m{l + 1 - m) + hlfl 
I m 



9{x, 



9i 



2ul 



m 



n h' 

ym"'l—m 



= 9ix)Y9nh{x)f, 



(39) 
(40) 



(41) 



(42) 



(43) 



where the equivalence signs in the above equations define new coefficients. In obtaining these, 
we have made repetive use of (^) . Note that the coefficient gf of Eq. ( |3^ ) is defined differently 
from (P). These notations are quite lengthy. However, they allow to obtain an algebraic 
equation very suitable for a symbolic computer program, generalization of Eq. ( |T7|) with (|l4|): 



N 



(N) 



ENi'^'^^flm) =0, 
|^AAr(AAr - 1) (^™°"^/°) Cfj^m^n) 



■m,j 



m=-N 
(TV) 



(44) 



fri) fj- 



2Xn (^'°''^- V°) f9]-„.-n + 2mA^ C^'^-'fn) fh] 



;'— m— n+1 



+ 



('ymiiiJ^O 



+ m 

l'ymin/^0^ /2i 1 



f'Vmin/0 



7°) h 



'j—m—n+l 



+ m{m - 1) (^-7°) (2/i]_^_„+2) - 2m (^™7°) ^]-m-n+i 

Vmin+Vmax 



+ 



ymin 



fc=0 



]-m 



(45) 



This algebraic equation has to be solved for each value of j. Apart from the possible general- 
ization using several basis functions (/ii(x), h2{x),- ■ ■) that we mentioned in the beginning of 
this Section, this equation is the most general algorithm we have found| to solve Eq. (||). 

In the next Section, we present several examples using both the generalized algorithm (|45|)- 
( p^ ) and the simplest one given by Eqs. (14) and (0). We use the latter when that the 



slight generalization is possible, consisting in taking, instead of the parametrization ( |24| ) of the potential, 
the ratio of two such paramatrizations. We do not consider this generalization here because we did not find an 
example needing it. This does not, however, exclude this possibility. 
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generalizations are not needed because it works faster. The examples are treated using a 
Mathematica implementation of Eqs. (|l|), (0), (|4|) and (||). 



4 Examples 

In this Section, we prove the convenience of our method for the determination of the eigenvalues 
and eigenfunctions of well-known Q.E.S. equations. We start with the genelarized sextic oscil- 
lator. We begin with a general set of parameters and show by fixing some of these parameters 
that the constraints between the different left parameters can be solved. In order to illustrate 
our method, a lot of details are given in the first treatment of the problem. Then, as we go on, 
only the main points are stressed. 

4.1 The generalized sextic oscillator 

The generalized sextic potential has the form 

V{x) = Vsb^x^ + {SYsab"^ + V2b^)x^ + {W^a^b + 2V2ab + Vib)x'^ 

^V,a' + V,u^ + V,a + V„ + ^^^;^ (46) 

with 6V3 > and ab > 0. It coincides with the sextic radial oscillator analyzed by Ushveridze 
iif 

Vs = f, V2 = 0,Vi = -At (^s + i + , Fo = 0, 

V-i = a(^s--^ (s-^^ ,y_2 = 0,a = 0,6 = 1, (47) 
while the Tkachuk potential is recovered for 

V, = 1, ^2 = ^1 = -3, V, = \a, y_i = \b, V.2 = \ab (48) 

in the context of two solutions. 

We see from (^) that V{x) can be written as Eq. (^) with the identification 

M = 3,/(x) =a + 6x2. (49) 

This, as well as f'{x) = 2bx, clearly suggests the basis function h{x) = x. From the rela- 
tions (|^-([l^) we get the nonvanishing coefficients: 

f^ = a,f^ = bJl = 2b,hl = l. (50) 

The unknowns of the problem are En, Xn, Cm \ 9i ■ They have to be found from Eq. (|T7|). 
Studying the case = 0, we can convince ourselves that the only nonvanishing coefficients of 
the weight function g{x) are (/ = 0, 1,2,3). Indeed, Eq. (17) reduces to 



^'JJ - Eq {a^5jfl + 2ab6j,2 + b'^Sj^i) = 0, j = 0, 1, 10 (51) 
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with — see Eq. (14) — 

4j = [-2ah\o-a^{glf + a^g\ + V-2 + ciV.i + a^VQ + a^Vi + a^V2 + a^V^5jfl 
+ [AabXogl " '^a^gWi + '^a'gl) 6j,i + [-462 ^ ^^^2^^ 
+ 4abXogl " 2a2(7,J<72 " a^glf - 2ah{glf + Sa^gi + 2ahg\ + + lahV^ 
+30^6^1 + Aa%V2 + 50^6^3] <^i,2 + (4a6Ao52 + 46^^050 
-2a^glgl - 2a^g\gl - Aabglg\ + Aahgl) 5j^^ + [iahX^gl + 462Ao5i 

- 2a2<7{<7i - a\g\f - Aabglg] - 2ab{g\f - b^g^f + Gabgl + b^gl + ftVo 
+3a6Vi + 6aH'^V2 + lOa^ftVa] 5^,4 + {'ib'^Xogl - "^a^dldl 

-^abglgl - Aabg\gl - 2b'glg\ + 2b^gl) 5j^5 + [^b^Ws - a\glf 

- 4abg\gl - 2ab{glf - 2b''glg\ - b\g\f + 2,b''gl + ftVi + 406^2 
+100^6^3] 6 + {-^abglgl - 2b'glgl - 2b'g\gl) 6jj 

+ [-2ab{glf - 2b'g\gl " b\glf + 6^2 + Sa^Vs] J,- g 

+ {-2b'gW3) ^9 + [-bHol? + b'Vs] ,5,, lo- (52) 

With vanishing g^,! = 5,... expansion coefficients, we would have get an extra term g^Sj^n, 
implying g^ = 0. Repeating this argument recursively, we see that only g^jl = 0, 1, 2, 3 survive. 
From the values j = 10, 9, 8, 7 we compute the weight function: 

g{x) = exp 




= eKpl^-^ayWsx'-\^V2x'-\V^,x'j, (53) 

while for j = 10 we have taken the positive root ^3 = ^/b-^Vs. We could have chosen the 
negative root. This would have lead to a nonphysical solution because it is not bounded when 
X ^ 00. We note however that unbounded solutions have their usefulness: they can be used 
to generate new potentials and the corresponding (bounded) eigenfunctions through Darboux 
transformations |14| . 

The index j = 6 fixes Aq- Together, the wavefunction ^0(2;) is determined to be 

ipoix) = Co{a + bx^) ^3 

X exp (-^a^/bVsx^ - ^^^^2^:2 - |V^^'^ , (54) 



the quantity cq being fixed by normalization. The next significant value j = 4 leads to the 
knowledge of Eq 



5 , 3 n 1 1 IV^ I V1V2 lb 

E, = ^„V3 + ^a'V, + -aV, + V„ - -o-^ + - - ^ -V,. (55) 
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Finally, taking care of j = 2 and j = 0, we are left with two constraints on the potential (46): 



IV-t IVlVi IV? 1 r-V? r- Vl 15, 1 r-V2 



VL 



64 Kf 8 



1 Vi 1 V1V2 1 2,. 

— a — ^ a a Vi 

16 Vi 4 ^3 8 

1 FiV? 1 

) = -a rf 

8 1-32 



16 



a^V2 



5 2^2^ 



15 



1 1 ^2^^ 1 2 1^1^2 

—a a — 5- a 

4 ^3 64 Kf 4 1/3 



32 Vz 64 

1 oFiFo 1 oV? 



(56) 



+ ^a^^ - — aVg 



16 - 



9 

64 



3 3 5 /- 21 1 31/2 5 /-K'' 

—a^V2 - -a^Vi - -aVb^^ ab H a^— + -aVb^ 

16 8 2 4 32 F3 8 ' 



2 
T 
•^3 



5,2 

4 



a'Vb- 



V2 



15 



(57) 



respectively. The different choices (47) (with n = 0) and ( ^q ) are compatible with these con- 
traints, which are limitations on the potential: There is at least one solution to the Schrodinger 
equation (||) with (^6|), this solution being given by Eq. (|5^, only if the different Vi^i = —2,..., 3 
coefficients of the potential satisfy the constraints. 

We proceed in a similar way for the = 1-case. Now, the set of equations (|T7| ) reads, for 
i = 0,...,11: 



.(1) 



(1) 



El {a^Sjfi + 2ab5j^2 + b^^jA) 



+ c 



(1) 



El {a^5j^i + 2abSj^3 + 6^5^,5) 



0. 



(58) 



The quantity h^^j is again given by Eq. (52) up to the replacement of Aq by Ai, while /i^^j is 

h'i^j = 2a^go6jfl + [-6abXi - a^g^ f + Sa^gl + V-2 + aV^i + aVo + aVi 

WV2 + 0^3] <5j, 1 + {4.ab\igl - '^a^oWi + 4a^52 + ^abgl) 6j,2 
+ [-Ab^Xj - 2b^Xi + AabXigl - 2a'glgl - a^{g\f - 2ab{glf + ba^gl + Qabg\ 
+V^ib + 2abVo + 2,a^bVi + Aa^bV2 + ba^V^] 3 + {^abXigl + Ab^Xig] 
-2a'^glgl - 2a'^g\g\ - Aabglg\ + Sabgl + 2b^gl) 5^-4 + [AabXigl + "^b^^wl 

- 2a^g\gl - a^{g\f - Aabglgl - 2ab{g\f - b^glf + 100553^ + 3b^gl + feVo 
+3a6Vi + 60^6^2 + lOa^^Vs] 5^-5 + {'ib^Xigl - '^a^dbl 

-Aabg^gl - Aabg\g\ - 2b^glg\ + Ab''g\) <5,- g + [462Ai53^ - a\glf 

- Aabg\gl - 2ab{g\f - 2b^glgl - b\glf + 5b^gl + b^ + 406^2 
+100^6^3] 5jj + {-AabgWs " '^b^aWs - ^b^abl) Sj,s 

+ [-2ab{glf - 2b'g\gl - b\glf + 6^V2 + 506^3] <5,- 9 
+ {-2b'glgl) Sj,io + [-bHgl f + 6^3] 

The analysis of the different values of j leads to the following eigenfunction 



ipiix) = Cl- 



X 



{a + bx'^y 



(59) 
(60) 
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corresponding to the energy 
iff the constraints on V{x) 



Ei = Eo + 2a V61^ (61) 



V^i{N = l) = V^i{N = 0)-aVb-^-a^^/b\^, (62) 

V^3 



V-2{N = 1) = V-2{N = 0) - Gab - aVb-^ + laVb^ 

V2 3 



"^3 



a'^VbVs (63) 



are satisfied. 

We have also computed the cases A'^ = 2, 3, 4. The results are too cumbersome to be written 
here. We have however observed that the constraints always change with N and that to ipN^x) 
corresponds the A'^th excited state (the identification is done looking at the number of nodes 
of the wavefunction) . Since the constraints change with N, this means that each obtained 
level corresponds to a different potential. This is at variance with the particular case a = 



(see Ushveridze Q). In this limit, we see from the comparison between Eqs. (54)-(p7|) and 
Eqs. (|6C|)-^(p^ that the cases = and = 1 coincide. 

We also see quite straightforwardly that the constraints (^)-(|5^) are compatible with the 
parametrization given in Eq. (47), provided n = 0. In fact, going to higher integer n, we can 



convince ourselves that there are n + 1 solutions compatible with the constraint on l^i. These 
solutions are defined on the half line ]0,oo[, due to the centrifugal barrier V-i/x"^, hence the 
name sextic radial oscillator. The higher n solutions are investigated using higher A^ values in 
our algorithm. 

In general, playing with the different parameters, we can generate a lot of different solutions: 

• In the general case above, we mentioned the fact that to one potential corresponds a 
single solution; 

• To the sextic radial oscillator corresponds a given number of solutions, depending on the 



value of the integer n in (47); 



Eliminating also the centrifugal barrier of the sextic radial oscillator, V^i = 0, we are left 
with a model defined on ] — oo,oo[, hence the name sextic oscillator. In that case, the 
exact solutions fall into parity classes. To see this statement, the parametrization ( |47| ) is 
replaced by 

V3 = t^V2 = 0, Vi = -4t (^n + I + , = 0, 

y_i = 0,F_2 = 0,a = 0,6 = 1, (64) 

where p denotes the parity. For a given n, one may get series of solutions: One with p = 
(even solutions), and one with p = 1 (odd solutions). Since p can take two values, the 
potential is then different for odd and even solutions: For a given potential, the solutions 
fall into parity classes, as said before; 
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Even with VLi and V-2 nonvanishing, we can have an exactly solvable potential. For 
example, the potential 

4 8 

^(^) = ^ + - n 2^2 (65) 

4 l + (l + x2)2 

is exactly solvable. This can be shown easily with our approach because we get no 
constraint. Theoretically, this is even more easily seen by noting that this potential can 
be constructed from the harmonic oscillator using a Darboux transformation [p 



The Kuliy-Tkachuk potential [R] obtained from the general case with the special parametriza- 



tion (48) is also interesting in the sense that it always generates a given number of solu- 
tions. This number is fixed from the start (for the Kuliy-Tkachuk potential, this number 
is three). 

In the following, we shall review the different possibilities we have just enumerated. 
4.1.1 The general case 

We have already discussed the general case given by Eq. ([46| ) for = and = 1. For higher 
A^, we get, before using the parametrization ( p7[ ) which still has to be shown to be consistent, 
the condition Aat = Aq — N/2. The equations for the constraints relating V-i and V-2 to the 
other coefficients, and the equations for the coefficients chi\ relating them to a single one to 
be fixed by normalization, for instance are too complicated to be solved analytically for 
arbitrary A^. However, we can find at least some of the solutions and see what they imply. Using 
our algorithm, one can see by direct substitution in the corresponding constraint equation that 

= compatible constraints, 

for which we also observe that E2N = Eq and E2N+1 = Ei. This however does not lead to 
new results: In fact, we just recover the case A^ = and A^ = 1, respectively. This is seen as 



follows, taking the case A^ = 2 as an example: Using the constraints (56) and (57), we can solve 

(2) (2) (2) (2) 

the equation for as a function of Cq . We get C2 — % b/a. Together with the relations 

Aat = Ao - 1 and f{x) = (a + bx^), we have fix)^''{c^o^ + c^q^x^) = /(x)^«4^Va- Then, the 
wavefunction 1^2 and ipo are equal, up to an unimportant constant. This analysis remains true 
to higher A^. This means that, when taking higher A^, we cannot take the same constraint 
equations for V^i and V-2 as the equations for A^ = or A^ = 1 to get new results. This 
implies that the constraints on V-i and V-2 are different for each value of A^, a result that we 
mentioned earlier. 

4.1.2 The sextic radial oscillator 

As we indicated above, this case corresponds to the parameter a = 0. From physical reason, we 
have to consider also V-2 = 0. This is in fact included in the constraint equations of the general 



case above, see Eqs. (57) and (|63|). To simplify the discussion, we choose, following Ushveridze, 
6 = 1,V2 = 0. As we have seen previously, the case A^ = and A^ = 1 are degenerate. We 
obtain 



+ (66) 



4VW 
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ipo{x) 
V-1 



3 4 
X 



Ei = 0, 
ipi{x) = coexp 



3 Vl 

2^ 4 4yTf 



(67) 
(68) 

(69) 



The last equation is taken from (56) or (|62| ) with the given a,6, V2 of this Section. It is clear 
that the parametrization ( |47| ) with n = satisfies it. 
Let us now look the case N = 2. We then have 



A2 

E2 



Ao-1, 



.(2) 



4e 



(2) 



(V^5)V4 



C2 exp 



3^4 W^2) 4 4VTi 



1 + 



ex 



= (1^3)3/4 



{Vl + 7^1^) (Fi + 9V^) 
4V3 



(70) 
(71) 

(72) 
(73) 



where e = ±1. We see that the last equation is satisfied by ( ^7| ) with n = 1. For this 
constraint, two solutions are possible: A ground state with e = 1 and a first excited state with 
e = — 1. Note that the argument under the square root is non- negative. This is seen easily 
with the parametrization ( [47| ) with n = 1. This can also be seen from the requirement that the 
wavefunctions must vanish at the origin, which implies —7/4 — T4/(4\/T^) > 0. This is enough 
to ensure a positive argument of the square root. 

The case = 3 leads to the same results as = 2. For A^ = 4, we get the new results 



A4 
E4 



An 



-4VV3 



.(4) 
.(4)' 



■04 (x) = Coexp 
with the three possibilities 

0,C4 



3 4 
X 



.(4) 



{Vl + uyVs) (Vl + 13V%} 



.(4) 2^3 



Vl + 8VV3 



.(4) 



i4)o.^/o y(-^l-9^g)^3 ^ ^(4) 



1^1 + 10^1 



-2c 



(4) 



2^3 



Vl + WVV3 



(74) 
(75) 

(76) 
(77) 

(78) 
(79) 



with e = ±1. This implies that we have three solutions corresponding to the constraint ([f7|), 
which is satisfied with the parametrization (|47|) with n = 2. These three solutions are the 
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ground state, obtained with ( [79[ ) and e = 1, the first excited state, obtained with (^), and the 
second excited state, obtained with ([7S|) and e = — 1. The arguments under the square root 
are non- negative. The argument follows the same lines as in the case n = 1. 
The case A*" = 5 is again identical to = 4. 

Using our algorithm, it is easy to see, and we have already mentioned it, that = Xq—N/2, 
or, because odd values of give identical results to the even values, X2n = Xo — n (see Eqs. (|70|) 
and (ff^)). This is true for arbitrary Vi and V3. From Eq. (|6^), a vanishing wavefunction at 
the origin implies — Vi/v^ > 4n + 3. This requires a negative Vi, as was also seen in the 
cases n = 0, 1,2 above. If we now choose Vi and V3 from (|47|), the condition is translated to 
s > 1/4, i.e., it becomes independant of n, and the exponent of is just s — 1/4, i.e., with the 
n-dependent parametrization of Vi, we have X2n = Xq = s — 1/4. 

From (p9D, ( |73| ) and ([77|), we can also find the condition that has to be satisfied by V-i: 

[V, + (4n + 3)^ [^1 + (4n + 5)^%] 

^-1 = Ws • 

With Vi and V3 from (47), we recover the corresponding V-i given by Ushveridze. 

Finally, the relation between the energy and the coeffcients Cm^ is E2n = — 4\/T^c^'^2/crr\ 
i.e., only the last two coefficients of the Taylor expansion of the wavefunction are required. 
This is clearly true when considering Eqs. (|6^), (|7l| ) and ([fsj). This is also true to higher order 
n. 

It is easy, although tedious, to go to higher n. For concision, the results are not written 
here. 

This Section has described the sextic radial oscillator. If we suppress the centrifugial barrier, 
the system is defined on the whole real line: The potential corresponds to the sextic oscillator. 
This is the subject of the next Section. 



4.1.3 The sextic oscillator 

As in the previous Section, we restrict the analysis to V2 = 0. We can solve the problem 
using the results of the sextic radial oscillator, implementing the constraint V-i = by taking 
s = 1/4 or s = 3/4, i.e., s = p/2 + l/A, with p, the parity, being equal to or 1. The 
parametrization ( |47| ) leads to (^4|). Compared to the radial case, we now deal with odd or even 
eigenfunctions. This is due to the fact that the potential is real and symmetric around the 
origin. The case s = 1/4 leads to a vanishing A^v, hence deals with the even eigenfunctions. It 
corresponds to p = 0. The case s = 3/4, corresponding to p = 1, implies Xn = 1/2. Starting 
with the sextic radial oscillator, defined on the positive half line, this leads to a prefactor 
/(j;)'^^ — Continuing the solution on the negative half line, we have finally generated the 
odd solutions. 

Since the wavefunctions correspond either to s = 1/4 or s = 3/4, i.e., to two different 
potentials, the solutions fall into two different classes: For a given sextic oscillator, it is not 
possible to generate simultaneously both odd and even solutions. 



4.1.4 An exactly solvable case 

Even with V-i and V-2 nonvanishing, the generalized sextic oscillator can be reduced to an 
exactly solvable model: Fixing the parameters in such a way to have V{x) given by Eq. (|65|), 
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all the solutions of the Schrodinger equation can be found. 

Using our algorithm (p7|) , we can see that, with f{x) = 1 + x^,h{x) = x, we get Atv = 
— l,g{x) = exp(— The general shape of the solution is 

^-(^) = rr^ ("t) ^"^^'^ ^^^^ 

with Cn a normalization coefficient and P„(x^) a polynomial of x'^ of order n. The first few 
values are: Pq = l,Pi = x{3 + x^),P2 = — 1 + + x^j-Ps = x{—5 + x"'). We see that a 
factor can be factored, with p = or p = 1, leading again, as expected, to odd and even 
eigenf unctions. Our approach also shows that the energies are related one to each other by 
En-\-i = En + 1, i.e., they are equally spaced like those of the harmonic oscillators. 

In fact, using an irreducible second order Darboux transformation, we can relate the po- 
tential (^) to the potential of the harmonic oscillator. This implies that the eigenvalues of 
the potential of this Section are the same (up to a shift) as the ones of the harmonic oscillator, 
excepted that we lose its two first excited levels. The Darboux transformation also relates the 
corresponding eigenfunctions: 

Po = 1, (82) 
P„+i = x(3 + x2)Fe„(x)-(l + x2)^^^, (83) 

with Hen{x) = 2~"/2jj„(x/\/2) where (x) is the Hermite polynomial of order n. 

It is trivial to check that the results coming from our approach satisfy this relation, and that 
our eigenvalues are the ones expected from the Darboux transformation. We note that since 
the Darboux transformation is constructed from the oscillator wavefunctions, the coefficients 
V-i, V^2 depend on the frequency lo. In the case we have investigated here, the values V-i = 
4, y_2 = —8 correspond to the particular choice Vi = 1/4, i.e., lo'^ = 1/4. 



4.1.5 The Kuliy-Tkachuk potential 

We now turn to another class of problems: Those for which the number of solutions is fixed from 
the start. As an example, we choose the Kuliy-Tkachuk potential Q, for which three solutions 
can be found. It is a particular case of the potential (|4^), but with V2 = V3 = 0. Before using 
the Kuliy-Tkachuk parameters given in Eq. (^), we first start with Vi,V(},V-i,V-2 arbitrary. 
Compared to the general case of the sextic oscillator, our algorithm allows us to see that only 
one coefficient Vi, say V-2, is constrained, the other ones being arbitrary. It also allows to see 
that the energies take simple values in terms of the parameters and A^: 

En = Vo + aVi + yWi [{2N + 1) + 4A7v] , A = 0, 1, 2, ... (84) 

However, as A^ is increased, the number of solution is growing because different values of the 
parameters are possible. For A^ = two solutions are possible: There are two possible values 
for V-2 in terms of a, b, V\. This leads to two possible values for Aq, and then to two energies 
and two wavefunctions. Note that they correspond to two different values of V-2, hence they 
do not refer to the same potential. The situation is similar for A^ = 1. For bigger A, we face 
the same difficulty as in the case of the potential (^6|): Having taken V2 = V3 = has not 
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simplified the solving of the problem, i.e., for general parameters, there is only one solution 
per given potential. On the other hand, fixing some of the parameters may lead to exactly 
solvable problems, such as in the previous Section, or to Q.E.S. systems with a fixed number 
of solutions : The K-T choice, with three solutions, corresponds to 



Vi = ^h, Vq={^- ^ 6, F_i = 2(3 - \/3)6, V.2 = (4\/3 - 6)6, a = 1. 



(85) 



With this parametrization, we find, for = 0: 



A„ = (80) 

Eo = 0, (87) 
/ a/3 a 

iI:q[x) = Co exp [ — —hx^ ) (l + hx^) i+v^ , 



while, for A'^ = 1, we have 



El = 3(2-\/3)6, (90) 
— —bx^\x{l + bx^)^+^. (91) 

The case N = 2 leads to the same solution as the previous cases = 0, 1, as well as to the 
new solution: 



2 



^(^-1), (92) 
E2 = 2(3-^/3)6, (93) 

Mx) = C2exp(^-^bx'^^ {l-bx^){l + bx'^)^^^~^\ (94) 

For higher A'^, we only generate the three cases above. As stated in the beginning of this Section, 
the Kuliy-Tkachuk potential corresponding to (^5[) allows to obtain exactly three wavefunctions 
and the corresponding eigenvalues which also satisfy the general relation (p^). 

Up to now, we have only considered the simplest variant (jl^) of our algorithm. We now 
discuss briefiy two simple cases were the general approach (|4^)-(^5|) is needed. These cases are 
the Lame equation and the screened Coulomb potential. 

4.2 The Lame potential 

The Lame equation is a Schrodinger equation with a periodic potential of the Elliptic Jacobic 
type: 

r ^2 ^ „ 1 

ipNix) = EN'ipNix), (95) 



+ k'^mi\]{mi^ + l)sn^(x, k) 



dx"^ 



18 



with k £ [0, 1]. This equation admits analytic solution for nonnegative integers m^r Q.E.S. 



extensions of the Lame equation can be found in the literature, e.g. [12|. In the following, we 
restrict ourselves to Eq. (|9^). This can be examined using our generalized algorithm (44)-(^). 
Several possibilities exist, showing the flexibility of our approach. For instance, we can choose: 



h{x) = sn(x, k) h'{x) = cn(x, A;)dn(x, k) = ib-y/[l — — k'^h'^(x)]. 

The lb sign implies to stay in a given interval, i.e. we choose the plus or minus sign in 
solution being extended appropriately on the whole line. Another choice corresponds to 



h{x) = exp(i(/)) h'{x) = ih{x)\ 1 + 



h{x) 



1 



h{x) 



(96) 
) the 

(97) 



with (j) the amplitude of x. A third choice, which is the one we shall take to illustrate our 
algorithm, is 



h{x) 



sn^(x, k) 



h\x) 



sn(x, k) 
cn(x, k) ' 

/l2(x) 
l + /l2(x)' 

+ h?{x){2 - k^) + h^{x){l - A;2), 



(98) 

(99) 
(100) 



to which corresponds f{x) = l/sn'^{x,k) = 1 + l/h^{x). Hence, we have Vmax. = 0, Fmin = 
2, V.2 = Vo = 0, V.i = k^mNirriN + 1), fj" = 6i^o + h = 8ifl, h = 6i,o + (2 - k^)5i^2 + (1 - 

For = 0, our algorithm gives immediately k = 0, or k = 1 and/or g{x) = 1. In fact, the 
two cases k = and k = 1 are included as special case of the full k problem, so that we can 
take immediately g{x) = 1. For a compatibility with the remaining equations, our program 
shows that we need either m = 0, Aq = 0, E = 0, or m = l,Ao = —1/2,E = 1 + /c^. The 
corresponding eigenfunctions are either a constant, or sn(x,A;), respectively. In fact, because 
Ao = —1/2, we see that the solution is il^{x) = [f{x)]~^^'^ = |sn(x,A:)|. However such a solution 
is not differentiable everywhere, and then is not a solution of the original Schrodinger equation. 
This implies that we need to consider the problem only on half a period of the function sn(x, k), 
thus allowing to remove the absolute operation, and to continue the solution on the whole real 
line, thus giving rise to ^{x) = sn(x, k). This need of studying the problem on a given interval 
was also present in the choice (|96[). The need was seen from the beginning because of the it 
sign at the level of h'{x). 

Putting the results together, we see that A^ = leads to 

(101) 
(102) 
(103) 
(104) 

or 



mo 


= 0, 


Ao 


= 0, 


Eq 


= 0, 


iI:q{x) 


= Co 



mo 



(105) 
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Ao 



1 

"2' 
1 + 



■ cosn{x,k). 

For = 1, several other new solutions are found: 



or 



mi 

Ai 
El 



nil 
Ai 
El 



1, 



1, 

cicn(x, k), 



-1, 

4 + k^, 

cicn(x, A;)sn(x, k), 



(106) 

(107) 
(108) 

(109) 
(110) 

(111) 
(112) 



(113) 
(114) 
(115) 
(116) 



or 



or 



mi 
Ai 
El 

ipo{x) 



2(l + k'^± \/l - A:2 + kA , 



ci 



(1 - k^)sn'^{x, k) + (k^^ Vl-k^ + k"^) cn^{x, k) 



mi 
Ai 
El 
V'i(x) 



3, 



3 

'2' 

5 + S/c^ ± 2\/4 - 7A:2 + Ak^ 



cisn(x, fc) sv?{x,k){k'^ 



I) 



2k' 



1 ^ V4 - 7A:2 + 4fc4) 



(117) 
(118) 

(119) 
(120) 



(121) 
(122) 

(123) 
(124) 



The solutions we got for = were also obtained. 

The analysis can of course be continued with higher values of A^. What we observe is the 
expected set [15| of eight different types of solutions modulated by polynomials in sn. 

4.3 The screened Coulomb potential 



The screened Coulomb potential ||16[ is also a nice example of the working of our method. It was 
shown in that reference that, under a given relation between the coefficients of the potential 
V{x) = F/x"^ + G/x + H/{x + z^), the solutions are polynomials. We shall verify this fact 



20 



using our formalism. We first take f{x) = (x + z'^) and h{x) = x with Vk = 5fc,o + <5fc,-i and 
= but for vifi = F5i_2 + GSi^^i and _i = H5ifl. We note that the coefficient F has to 

be -F > —1/4. This comes from the physical requirement of having, with the parametrization 

F = L[L + 1) + Fi and L a nonvanishing integer (angular momentum), Fi > — (L + 1/2)^. 
Starting our algorithm with = 0, we obtain straightforwardly that one of the coefficients 

of the potential must be expressed in terms of the others, e.g, G: 



G 



4F + (1 + VI + 4F) (4 + H2 



2z2 



Then, the solution can be expressed as 



Ao 

Eq 



1, 



4F + 2 (1 + VI + 4F) (1 + Hz^) + {Hz^) 



2\2 



4^4 



exp 



^ (l + VTTlF + Hz')] ^(i+^T+37)/2 + 



(125) 

(126) 
(127) 
(128) 



with a positive argument of the square root from the condition F > 1/4. Another solution is 
possible but we have rejected it because it is unbounded at the origin. We also observe the 
following relation between the energy, and the different coefficients of the potential: 



(G + HY 
-En 



3 + V1 + 4F 



(129) 



In fact, the vanishing of the wavefunction at infinity requires Hz'^ < 1 + Vl + 4-F, see ( |128 ) 
Because of this, it is easy to verify that (G + H) < 0. Hence, Eq. (129) can be replaced by 



(G + H) 
V—Eo 



3 + V1 + 4F , 



(130) 



a relation which was also obtained in |16|, with the index M = 1 in that reference. 

The analysis can of course be done for higher N. The case = 1 is however already 



highly involved and will not be presented here. We just quote that, as in [16|, the number 3 in 
Eq. (pO) has to be replaced by 2M + 1, with M > 1. 



5 A new Q.E.S. potential 

Using a Darboux transformation on the following sextic radial oscillator 

Vix) = \x^-5x' + ^, 

we can show |14] that the wave function 

exp (-x^/S) x3/2 (6 + x^) 



V'(x) 



c8 + 4x4 + 20 



(131) 



(132) 
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is a solution of vanishing energy of the Schrodinger equation with the potential 



N 1 6 2 3 16x2 (x^- 6) 2048x6 , ^ 

y x = -x^ + x2 + — + — ^— L ^. 133) 

^ ^ 4 4x2 x8 + 4x4 + 20 (x8 + 4x4 + 20)^ 

We show now this proposition using our algorithm, noting first that the potential obeys the 
parametrization /(x) = x^ + 4x^ + 20 and h{x) = x. Moreover, Vk = 5fc,o + ^k-i + 5k-2 and 
vi^k = but for vifi = 3(5/ _2/4 + (5/,2 + 5i,6/4, u/ -i = -96(5/,2 + 165i,6 and _2 = -2048(5/,6. 

The cases = and = 1 do not lead to a solution. For N = 2, all the constraints are 
fulfilled provided that g} = l/8Si^4,gl^ = 1'^^ together with cL^) = cf^ = cf^ = 0,cL^^ = 
1/4(17 - 25/„g)4^^ = 64^\ and E2 = 0, A2 = l/8(5/„g - 9/2) = -1. The wavefunction is then 



^^{x) = cX' exp ( -— \ \ ^^x']{xf' ^ 4x^ + 20) , (134) 



which is nothing else than ( |132D . 

6 Conclusions 

We have proposed a general method for determining the analytic eigenfunctions (and associ- 
ated eigenvalues) of a given (quasi) exactly solvable Schrodinger Hamiltonian. This method 
consists in a general algorithm replacing the differential Schrodinger equation by a finite set 
of algebraic equations, these ones being treated through a Mathematica implementation. Our 
general algorithm might be considered as relatively heavy from a purely technical point of view 
but it actually works straightforwardly. Moreover, for the majority of analytically solvable 
Schrodinger equations, it can be replaced by a simpler version as presented in Section In all 
cases, it gives rise in a very systematic way to the solutions of the Schrodinger equations and 
also precises the conditions under which these equations are solvable. This systematization is 
particularly useful at the level of quasi exactly solvable Schrodinger equations for which one 
knew before our approach that some of them are relevant of s/(2, K) and the others not. These 
two different points of view have been unified through our formalism. We also assert that new 
quasi exactly solvable potentials (see Section ^) can be handled within our algorithm due to 
the fact that the group theoretical approach is not needed anymore. 

We are also convinced that it is possible to extend our algorithm to the relativistic context, 
studying more particularly the quasi exactly solvable Dirac equations. Only a few studies, 
see for instance Ref. |17|, deal with these equations because, by opposition to some of the 



Schrodinger ones, they are not in general subtended by a specific Lie (super) algebra. Our 
approach avoiding this group theoretical aspect should be of interest in the analysis of such 
Dirac equations. We plan to come back on these developments in the near future. 
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